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SPECTRUM OF TURBULENCE IN A CONTRACTING STREAM! 
| By H. S. RxsNER and M. TUCKER 


SUMMARY 


The spectrum concept is employed to study the selective effect 
of a stream contraction on the longitudinal and lateral turbulent 
velocity fluctuations of the stream. By a consideration of the 
effect of the stream contraction on a single plane sinusoidal 
disturbance wave, mathematically not dissimilar to a triply 
periodic disturbance treated by G. I. Taylor, the effect on the 
spectrum tensor of the turbulence and hence on the correlation 
tensor is determined. Lack of interference between waves 
follows from the postulation of a low level of turbulence; this 
and the assumption of an inviscid fluid imply neglect of decay 
effects. The compressibility of the main stream is taken into 
account, but the density fluctuations associated unth the turbulence 
are assumed to be negligible; this would be the case 1f the turbu- 
lence originated from wakes and boundary layers in the very 
low-speed portion of the flow. For an axisymmetric contraction 
and a particular isotropic initial turbulence, some explicit 
results are obtained. The one-dimensional longitudinal 
spectrum is found to be.distorted (as well as reduced in ampli- 
tude) with its peak shifted well to the right of the initial position 
above the zero of the wave-number scale. The selective effect of the 
contraction on the mean square longitudinal and lateral com- 
ponents of turbulent velocity 4s found to be gwen uniquely when 
the initial turbulence is isotropic, regardless of the details of 
the spectrum. If the initial spectrum 18 anisotropic, as, for 
instance, that produced by a damping screen, then the selective 
effect is altered. 

In a crude extension, decay effects outside the scope of the 
theory are allowed for in first approximation. With this 
extension, a comparison with experiment is made of the selective 
effect on turbulent intensity where the estimated decay effects 
are comparable with the contraction effects. The agreement is 
good for the longitudinal component, very poor for the lateral 
component, the experimental data themselves being in conflict. 


INTRODUCTION 
The generation of a wind-tunnel flow is always accom- 


panied by a certain amount of turbulence; this is one respect ' 


in which the flow fails to simulate free-flight conditions. 
Measurements in the tunnel, particularly those sensitive to 
boundary-layer behavior, are known to be affected by this 
turbulence. Accordingly, the tunnel designer attempts to 
reduce the intensity to the lowest practicable level. The 
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use of honeycombs and damping screens in a large low-speed 
section (settling chamber) followed by a sharp contraction 
to the much-higher-speed working section is known to be 
effective. The honeycombs and screens located in a low- 
speed section reduce the absolute level of the turbulence 
with little drag penalty; then the relative level is greatly 
reduced by the large gain in tunnel speed through the con- 
traction, aside from any effect a the contraction on the 
absolute level. 

Once the characteristics of ihonepoombe and screens are 
known, the further quantitative estimate of the reduction in 
turbulence involves a knowledge of the effect of the tunnel 
contraction? on the turbulence. It is known that the 
longitudinal component of the turbulence is greatly reduced 
(in absolute value) by the contraction; the behavior of the 
lateral component appears, on the other hand, to vary from 
no change to a substantial increase. Prandtl (ref. 1) ob- 
tained a quantitative estimate of the first effect by consider- 
ing the conservation of energy for a perturbed longitudinal’ 
filament: if the initial stream speed is U, the filament speed 
U-+u, with 4«U, and the final stream speed is LU, then 
the final filament speed must be /U-+-l-‘u; that is, the con- 
traction reduces the longitudinal perturbation velocity u, by 
the factor (7. For the lateral effect, Prandtl applied con- 
servation of momentum to a small rotating cylinder of the 
fluid, with its axis cross stream, as the fluid traversed 
the tunnel contraction. He concluded that the lateral 
perturbation velocity v is increased by a factor +l. 

Prandtl’s considerations on the effect of a stream con- 
traction were limited, as has been noted, to particular 
idealized ‘‘eddies.” G. I. Taylor (ref. 2) later attempted 
more realism by treating a mathematically defined model of 
turbulence which amounted to vortices in parallelepiped 
partitions arranged in a regular three-dimensional array. 
The changes in vorticity on traversing the contraction were 
determined from a theorem based on conservation of circu- 
lation for an inviscid fluid; the corresponding altered turbu- 
lent velocity pattern was then calculated. The final result 
of the analysis consisted in expressions for the root-mean- 
square longitudinal and lateral turbulent velocity compo- 
nents wu’ and v downstream of the contraction expressed as 
ratios of the corresponding values upstream. 
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The initial condition of isotropic turbulence (mean values 
unaffected by rotation or reflection of axes) was approxi- 
mated by specifying the vortex partitions to be cubical. For 
this case the reduction in the longitudinal component u’ was 
found to vary more nearly like 1.5 /7! than the value /-! sug- 
gested by Prandtl. No explicit result was found for the 
variation of the lateral component, however: the calculations 
contained a free parameter. 

Taylor's results for the longitudinal component agreed 
fairly well with the experimental data then available, but it 
is now considered that the measurements were made too close 
behind the screens for the screen-produced turbulence to 
have been isotropic. On theoretical grounds, the objection 
to Taylor's theory is threefold: first, the decay processes of 
turbulent mixing and viscous dissipation, which result in a 
reduction of the mean intensity with axial distance in the 
wind tunnel, are neglected; second, the assumed model of 
turbulence fails to exhibit the spatial and temporal random- 
ness of actual turbulence; third, no choice of the parameters 
in Taylor's model corresponds to isotropy. In a sense all 
three objections apply likewise to Prandtl’s results: no model 
was employed in his considerations, and hence no distinctions 
between the effects of isotropy and anisotropy were made. 

The second and third objections can be removed by work- 
ing, not with a model of turbulence, but instead with a 
Fourier integral representation of & random turbulent field. 
The integral can be interpreted as a superposition of plane 
transverse sinusoidal waves of all wave lengths and with 
apparently random orientations and phases. This aggregate 
of plane waves constitutes the (three-dimensional) spectrum 
of the turbulence. Only the statistical aspects of this 
spectrum will be known, not, for example, the detailed phase 
relations. Mean-square velocity components "may be ob- 
tained by an integral of certain spectrum functions in which 
the phase relations are suppressed; these functions are in- 
cluded in the ‘spectrum tensor’ (ref. 3). 

Taylor’s concepts may be applied to find the effect of a 
stream contraction on a single plane wave. The effect 
under the assumptions is linear; therefore, the superposition 
implied by the Fourier integral may be employed to obtain 
the contraction effect on a field of turbulence. In particular, 
if the initial spectrum tensor is known, the final spectrum ten- 
sor is determined; the initial and final mean-square velocity 
components then result from quadratures. Indeed, from 
` the same information, changes due to the stream contraction 
in correlations of velocity at different points may be calcu- 
lated: use is made of the fact that the correlation tensor 
is an inverse Fourier transform of the spectrum tensor 
(ref. 3). 

Accordingly, the injection of the spectrum point of view 
into Taylor’s original concept of the contraction effect 
makes possible a more realistic calculation of the changes 
in mean-square velocity components. In addition, it 
provides much more detailed information concerning 
changes in the statistical properties of the turbulence; that 


is, it provides the changes in the spectrum tensor and in the 
correlation tensor. l 

The ideas just outlined are developed in the present paper. 
The first section is devoted to an account of turbulent 
spectrum analysis in a form specially adapted to the analysis 
of the contraction effect. In this account, which is a general- 
ization of a development in reference 4, the role of the spec- 
trum tensor is subordinated to that of the individual Fourier 
components (plane waves) in contradistinction to the cus- 
tomary treatment. This approach has perhaps an auxiliary 
merit in providing some better physical insight into the 
significance of the spectrum tensor. 

Next the effect of a stream contraction on a single plane 
wave is calculated by an application of Taylor’s concepts. 
The treatment is slightly more general in that compressibility 
of the main stream is allowed for. The density fluctuations 
associated with the turbulence are assumed to be negligible; 
this would be the case if the turbulence originated entirely 
from boundary layers and wakes in the very-low-speed portion 
of the flow. Following Taylor, the problem is linearized by 
postulating a sufficiently weak turbulence so that the self- 
distortion of the turbulent eddies is small compared with the 
distortion imposed by the contraction of the main stream; 
this together with the assumption of an inviscid fluid m pitos 
neglect of the decay of the turbulence. 

In succeeding sections the spectrum and correlation tensors 
downstream of the contraction are expressed in terms of 
the corresponding initial tensors. For the special case of 
an axisymmetric contraction and isotropic initial turbulenco, 
the ratios of the root-mean-square longitudinal and lateral 


velocity fluctuations downstream and upstream are obtained 


explicitly in terms of the parameters defining the contraction. 
For a particular subcase where the initial isotropic spectrum 
tensor is specified, the corresponding ‘one-dimensional’ 
spectrums (as would be recorded by stationary hot-wire 
probes) upstream and downstream of the contraction are 
calculated; the specification is such that the upstream one- 
dimensional spectrum corresponds to experiment (ref. 5, 
p. 35) in a number of cases of isotropic turbulence. 

Most of the calculated contraction effects are amenable to 
experimental checks either directly or indirectly. The 
available experimental data, however, are limited to the 
changes in the root-mean-square velocities. A comparison 
with these experimental data is given with an estimated 
allowance for decay effects’ outside the scope of the theory. 
Design curves of the changes in the root-mean-square velocity 
components neglecting decay are included for engineering 


urposes. 
E SPECTRUM ANALYSIS 


Representation of turbulence by superposition of plane 


‘sinusoidal waves.—Suppose di, 4s, gs represent the com- 


ponents of velocity in a turbulent field; that is qı, qa, and qs 
vary in an apparently random manner in space and timo, 
and the mean values 4,—3;—45—0. Subject to certain 
conditions, a snapshot of this field at any instant can be 
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represented as a set of three-dimensional Fourier integrals 
Qa (21, t3, :3- | | f e. (Ei, ka, ka) ef Circ hr bord dkidkadk, 


where a=1, 2, or 3 and the significance of k, ka and ks 
will be brought out later. A continuous representation of 
the turbulent field is obtained by allowing the Qa to vary 
with time. 

It will be convenient to abbreviate the Fourier integral to 


qa(x)= | f [Xo (18) - 


and to introduce the companion equation 


Q.0)—-8« | f [2.747 (9) "by 


where k—4, ka, k3; 2=2%, x, 2%; dr (x)=—dz,dx,dx;,. The 
second equation allows, in principle at least, the coefficients 
Qa(k) to be calculated. Mathematically, q«(z) and Q.(%) are 
termed three-dimensional Fourier transforms of each other; 
use will be made of this relation later. 

The velocity components qa are connected by the condition 
of continuity. In many cases of practical interest, these 
turbulent velocities originate from boundary layers and the 
wakes of obstacles in flows of low subsonic speed, so that 
associated density fluctuations may be ignored; this 1s still 
permissible when the turbulence so produced is transported 
by a high-speed stream. ‘Thus the incompressible form of 
the continuity equation may be used and the result is 


Qe + Qakta+ Qk; =0 i 


This relation may be written more compactly as 
2 Qaka=0 (2) 


Physical interpretation.—The amplitude components Q. 
are complex in general. According to equation (1b), then, 
the requirement tbat the velocity components ga be real 
implies that Q.(—£) is the complex conjugate of @.(k). If 
corresponding terms for k and — k in equation (18) are paired, 
their sum is thus equal to the real quantity 


2(Re Qa) cos (k-z) —2(Im Qa) sin (k-z) (3) 


The imaginary parts cancel in the pairing, which implies that 
they contribute nothing to the integral. Expression (3) 
represents a pair of plane standing waves, a cosine wave and a 
sine wave, with normals in the direction k=(k, ke, ks), 
where z= (zı, %2, 43) is the radius vector to any point. The 
vector k is termed the wave-number vector and its magnitude 
k simply the wave number; the corresponding wave length is 
2r divided by the wave number. Since k is perpendicular to 
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the wave front, it is sometimes referred to herein as the ‘wave 
normal.’ f 

The continuity condition, equation (2), states that both 
the real part (Re Qa) and the imaginary part (Im Qa) of the 
amplitude vector Q— Q,.— (Qi, Qe, Q4) are perpendicular to 
the wave normal k; that is, both waves of expression (3) are 
transverse. For each wave any one of the parallel planes 
containing beth the local velocity vector g and the wave- 
number vector & (which is perpendicular to Q and hence to q) 
is called the plane of polarization. The cosine wave (real 
part) and sine wave (imaginary part) may be polarized in 
different planes in general; the necessary and sufficient con- 
dition that they be polarized in the same plane is 


ReQi ReQ: ReQs (4) 
Im QO, Im Qi Im Qs 


Equations (1) are now seen to represent a superposition of 
plane sinusoidal waves (Fourier components) with all orienta- 
tions of the wave fronts (all directions of the wave normal E) 
and all wave lengths (all wave numbers k). Each wave is 
transverse, and all planes of polarization. are permitted. 
For each value of k there exist a cosine wave and a sine 
wave; their respective amplitudes and planes of polarization 
are different in general. The complex amplitude components 
Q.(k) express, in their real and imaginary parts, how the 
respective amplitudes and planes of polarization vary with 
the wave-front orientation and the wave number. 

Mean values: correlation tensor.—Consider the spatial ? 
mean value of the product of the velocity component qa at x 
and the velocity component qg’ at z'/—z-L-r as x varies but 


- the separation r of the two points remains fixed during the 


averaging process; this mean value is called a velocity cor- 
relation and is given the symbol E,s(r). There are nine such 
correlations, corresponding to a=1, 2, 3, and B=1, 2, 3. 
The form Z,s(r) has been shown to transform like a second- 
order tensor and has been designated (sometimes divided by 
q7) as the ‘correlation tensor’ (ref. 6). 

Evaluation of correlation tensor in terms of spectral 
quantities.— The basis of the idea that the correlation tensor 
might be expressed in terms of individual-wave parameters 
is drawn from reference 4. In that paper mean-square values 
of gi, d», q in terms of such parameters are discussed; these 
constitute the diagonal terms of Z.,s(r) evaluated at r=0. 
The following derivation of R.s(r) amounts to a generaliza- 
tion of the derivation of #,,(0) given in that reference. 

Assume that the turbulence is confined to a large, but not 
infinite, parallelepiped of edges 2D,, 2D», 2D; and vanishes 
everywhere outside. The space average Ras(r) is to be 

3 If the statistical properties of the turbulence are independent of position (homogencous 
turbulence) and time-independent, an average at a given time over all space equals an average 
at a given point (or pair of correlated points) over ali time; a proof is given in appendix B. 
If the statistical properties vary slowly with time, the space average will still approximate a 
time average over an interval just long enough to smooth out the fluctuations. 

4 The space integral of the square of any component velocity is thereby bounded, which isa 
prerequisite for the existence of a Fourler integral, equations (1); in other words, equation 
(1b) shows that Qa must depend on the volume 8D; D:D, of the region within which ga differs 


from zero: if this volume is infinite Qa Is infinite, and the Fourier integral, equation (18). does. 
not exist. 1 
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taken for all points z in the interior, the separation r remaining constant and very small compared with 2, ; 3. 


^ A jJ ff [eionanennn es Án 


MEE Ds 
Rar) q«(z) qe(x Tsp ool. 


or, interchanging the order of i EO 


The average is 


Put Ay=kh, +k’ ; P } and m -h E. 
kı ark Tat esr; 


Then recalling that k 


=(k,, k ks), dr(k)=dk,dkadks, etc., and ` k’r= 


Bo D E; f f f Qs les) í f f f QU, — kı, Fs — ka Ea — ka) 


pH UC -tratr Gs rg SD Ki Dy sin KaD; sin eS dE, AES d Ks, 5 dkıdkadkı 


The inner integral approaches a limit as the parallelepiped edges 2D,, 2D,, 2D; become indefinitely latge. 


limiting case tbe correlation tensor is then 


ED ONE CM 
For thia 


Rea(t)=lim f f f DET. Q. esee Qk ke, — ege td dd eg lim f f f ST QEG Erdr(E) (6) 


where Qs* (k) I= Qs(—£)] is the complex conjugate of Qa(k) and 7 is the volume 8D,D4D,; of the parallelepiped. 


Identification of lim = — T 9.Q, with spectrum tensor.— 


If the field of duties is homogeneous, the correlations 
should be independent of the volume 7 averaged over when this 


3 
volume is sufficiently large; thus lim ŠT. Q«.Qs* should exist. 
From the mathematical standpoint, equation (5) shows that 
Ras(t) is the Fourier transform of lim 2*- Q,(k)Qs*(E), and 


conversely; since it is known that the integral of |Rag(r)| 
over all r is finite, the inverse Fourier transform relation 


3 ` 
ensures that lim on Q«.Qs* exists. oa since the Fourier 


T QuE) Qs*(k) may 
be identified with-the form T'as(k) defined = Batchelor (ref. 3) 
as the Fourier transform of Aas(r). The form Tas(k) is 

known as the spectrum tensor. The Fourier transform rela- 
tions connecting the spectrum tensor and the correlation 
tensor are summarized as 


transform of à function is unique, lim == 


TÉ) 5s zi J [Baoi 0 


Ral) = j f [Ttt 00 (6b) 


By use of the Fourier transform relation, Batchelor demon- 
strated that Tas is a second-order tensor and obtained a 
number of its properties. Thus, for example, Tag is complex, 
in general, with Tga=Tas*, and the diagonal elements Tea 
are real; also, Tas(—k)=Tga(k). It is of interest to observe 


that these same properties result immediately from the 
identification of lim oe Q.Ve* with Tas. Thus Q.Qs* is 


complex, in general; QsQ;* equals [QaQs*]*; and Q&Q.* is, 
of course, real. Furthermore, since Q.(—k)=Q.*(k), 
Qa(—k)Qe*(—k) equals Q.*(E)Qs(k); hence Tas(—k)= 
Pga (£) . 

The distinction between cases where I'g is real and cases 
where it is complex may be given & physical interpretation. 
The product (,95*, and correspondingly Tag, is seen to be 
real when the condition equation (4) is satisfied. This 
implies that the cosine wave and sine wave associated with 
wave number k are polarized in the same plane. The alter- 
native condition where (Q4,05*, and hence Tap, is complex 
implies polarization of corresponding cosine and sine waves 
in different planes. The velocity pattern of such a pair of 
waves is quite interesting: successive velocity vectors along 
a line in the direction of the wave normal & turn progressively 
about this line in spiral fashion; the tips of the vectors trace 
out a helical curve on a cylinder of oval cross section. 

Energy spectral density.—Each of the diagonal elements 
Iu, Da, and Pa of the spectrum tensor Tes may be inter- 
preted as an energy spectral density. Thus, according to. 
equation (6b) 


Ry (0) =w = f f fracas (k) 


f 
Therefore, the differential 5 luüdr(É) represents the contribu- 


tion to the kinetic energy component 5 wj per unit mass 


made by waves with wave number within the range dr(k). 
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One-dimensional spectrum.—The elements of the three- 
dimensional spectrum tensor are not directly measurable; 
they may be obtained by taking the Fourier transform of the 
measured correlation tensor. A hot-wire probe placed in 
the moving stream will, however, develop a fluctuating out- 
put voltage whose (one-dimensional) frequency spectrum 
(ref. 7) is related to a diagonal element of the three-dimen- 
sional spectrum tensor. Thus by equation (6b) the contri- 
bution to the mean-square velocity component qa (=Raal0)) 
from all waves with wave-number components in the direc- 
tion of the z-axis between |k,| and |k,|+-|dk,| is 


FiAkidk, =2 (J fritar) dk, | (7) 


the factor of 2 accounting for suppression of negative values 
of kı. The function Falki) is the one-dimensional spectrum 


corresponding to the velocity component q.; the values 


a=1, 2, 3 correspond respectively to the longitudinal and 
two lateral spectrums. The particular spectrum obtained 
depends on the arrangement of the hot-wire probe elements. 
EFFECT OF STREAM CONTRACTION f 

Consider now that the turbulent velocity pattern qi, qs, q 
is carried along by an inviscid general stream with velocity 
U(a,) in the x-direction. Consider also that qi, q2, q3 are 
so small that their effect on the streamlines may be neglected 
as the flow traverses a wind-tunnel contraction. The con- 
traction will, however, distort the shape of fluid elements. 
(See fig. 1.) The vorticity distribution will be forced to 
alter accordingly to conserve the circulation about each 
element. The net result will be an altered pattern of turbu- 
lence. Each plane wave (Fourier component) Qi, Qo, QaeE z 
will, in fact, be altered independently under the linearizing 


assumption to be made; the. over-all effect on gı, qs, qs can . 


be obtained by the summation expressed by equation (la). 
It thus suffices to consider the effect of the contraction on 
a single representative plane wave. 


EFFECT OF CONTRACTION ON REPRESENTATIVE PLANE WAVE 


Velocity and vorticity at upstream station.— Designate by 
A a reference station upstream of the contraction and by Ba 
reference station downstream of the contraction. (See 
fig, 1.) Let & typical Fourier component (plane wave) of 
the turbulent field g«(a—1, 2, 3) at station A be represented 
at time t=0 by 


at= Qat €? (8) 


This wave, equation (8), is supposed to be carried along by 
the main stream with velocity U. 


The vorticity o4 is obtained from the curl of equation (8) 
08 


G.A i Leg ig Qu tto (9) - 


0, if any pair of subscripts are equal 
where S , if aBy are in cyclic order 
—1, if «By are in anticyclic order ^ 


f 
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(b) Stream tube geometry. 
FIGURE 1.—Schematic representation of flow contraction parameters. 


Distortion of fluid element in passage through contrac- 
tion.—Suppose the contraction is such that the stream 
velocity U is increased by a factor J, between stations A and . 
B while the breadth and the height of the tunnel are reduced 
by factors 4 and l, respectively. (See fig. 1 (a).) In 
traveling from A to B, an initially cubical element of fluid 
of edge D will be distorted into a parallelepiped of edges 
LD, LD, GD (see fig. 1 (b)); a particle in the element original- 
ly (t=0) a vector distance z from a corner particle will finally 
(£—1) be found a distance £ from the corner particle, where £ 
is related to z by ` 


£,=1,2, 
Eb, (10) 
f= hts 


In this argument the modification of the streamlines due to 
the turbulent velocity fluctuations has been neglected. 
This implies that the relative displacement of two adjacent 
particles due to the superposed turbulent motion is small 
compared with the displacement due to the tunnel contrac- 
tion. This key assumption, due to Taylor (ref. 2), linear- 
izes and vastly simplifies the problem. The limitations 
imposed by the assumption are discussed later under DECAY 
CONSIDERATIONS. 

The velocity ratio l, and the lateral and vertical contrac- 
tion ratios 1, and J; bre related by the continuity condition 


al4lsl3— 1 


i 
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where c is the ratio of stream densities at stations B and A; 
the density is considered uniform at each station in accord- 
ance with the initial assumption of negligible turbulent 
density fluctuations. 

Vorticity at downstream station.—The vorticity is carried 
along by the flow, the fluid elements undergoing the distor- 
tion pictured in figure 1 (b), to the approximation used. 
During the motion the strength changes in such a way as to 


maintain the constancy of circulation of the fluid elements.’ 
The changes are expressed by the equations for the transport 


of vorticity in the Lagrangian form, due to RI (ref. 8, 
p. 205). 

49£« 
dg Og 


where o is the density ratio between stations B and A, and 
the derivatives d£,/Oxg express the effect of the fluid distor- 
tion. Evaluation by means of the distortion equations (10) 
yields simply 

Ge? = laa" 
or, in expanded form i 


Q j= g Lo i 


wP — clo, 
@3°= clwt 


'These equations relating downstream and upstream vorticity 
embody the entire dynamics of the contraction effect. The 
equations are not limited to the plane sinusoidal waves dis- 
cussed earlier, but apply to any (weak) vorticity distribution 
whatsoever. 

The above derivation of the vorticity changes 1s substan- 
tially in the form given originally by G. I. Taylor (ref. 2) for 
the case c=1 (incompressible flow). In order to assess the 
influence of the simplifying assumptions a more general 
derivation based on the Navier-Stokes equations is given in 
the section entitled DECAY CONSIDERATIONS. 

By virtue of equation (9) as applied to equation (lia), 
the vorticity at station B is obtained explicitly as 


a iol > | Eagrkg yetz 
Y 


(118) 


(11b) 


where, it will be remembered, x is the radius vector to a 
fluid particle at time ¿=0 when the fluid element is at 
station A in the moving coordinate system of figure 1 (b). 
The corresponding vector to the particle at time t=t, when 
the fluid element is at station B, is £ in that figure. When 
equations (10) are used to express z in terms of £, the ex- 
ponential k-z becomes, in expanded form, g 


k-z 


ted inviscid fluid. The modifica Mon proaneea i by 
case of a gas, becomes appreciable for the smaller ed 
a criterion for neglect of viscous Ties. P 


3 This statement is exact for the 


her wave numbers; for 
»k1««ldU]dr]. (See DECAY CONSIDERATIONS, deuation (43) > 
è These equations refer to axes moving with some fluid particle rather than axes fixed as in 
reference 8; the form of the equations is unaffected. 
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The right-hand side may be expressed as x -£, where 


E, ka ks 
db 


defines a new wave-number vector. 

Velocity at downstream station, general case.—Tho 
velocity distribution whose curl in the £j, &, & system is 
given by equation (11b) and which satisfies continuity is 
found to be expressible in the form 


GaP =Q Fett 


QP =F ( 0.4— a) 


lg kt 
where x is the magnitude of the wave-number vector x. 
This result is the general solution for the contraction effect 
on a single plane wave. 
Equations (12) and (13) admit of a simple (but not obvious) 
geometrical interpretation: traversal of the stream contraction 
alters the initial plane wave, equation (8), so that tts wave- 


with 
(13) 


number vector k—(k,, ks, kz) is transformed into ==- [^ 
1 


and its amplitude vector (Q4, Q4, Q4) is transformed into the 


projection of (Q,4/l,, Ql, NA on a plane normal to the : 


new wave-number vector x. 

Velocity at downstream station, axisymmetric contrac- 
tion.—In case the stream contraction is axisymmetric,’ a con- 
siderable simplification results. The condition for axisym- 
metry ,=/;, with use of the continuity equation (2), reduces 
equations (13) to 


6 P= _ Qt bà EST ES 
d ek? Mr SEES 
A Q^ kık(1— e6). 
e | ate ek tbh Pky 


um 1 Qik, E. = e) 
B ls | Qs ERES 


(14) 


where ¢=],7/l,7. The considerably greater complexity of 
equation (13) is perhaps obscured by the purposely expanded 
form of equations (14). 

Tf the initial wave normal k is perpendicular to the (longi- 
tudinal) z,-axis, the component £, vanishes and equations 
(14) reduce further to 


= QiA/L 
Q7 = QsAIl; (15) 
Q= Qs4/ ls 


The same equations result when QA may be neglected in 
comparison with Q;4 and Q4, that is, when the amplitude 
vector is substantially normal to the z,-axis. Equations (15) 
state that an axisymmetric contraction defined by lı, la alters 


7 A contraction such that all cross sections of the tunnel aro similar, whence 13315, is termed 
arisymimetric; the sections need not be circular. 


(12) 


~ 
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these waves by a factor of 1/l, in the longitudinal velocity comi- 
ponent and a factor of ifl in the lateral velocity components. 
These equations apply only to particular types of wave; yet when 
the contraction effect 1s later integrated over the random aggrega- 
tion of waves representing rsotropic turbulence, the over-all 
resulis are found not to differ greatly from the simple factors 
1/I, and 1/le, respectively. 

The same factors were obtained by Prandtl (ref. 1) for 
other special disturbances: the factor 1/l, from energy con- 
siderations for purely longitudinal disturbance velocities, and 
the factor 1//; from momentum considerations for a rotating 
cylindrical element of fluid with axis normal to the stream. 


EFFECT OF CONTRACTION ON SPECTRUM AND CORRELATION TENSORS 


Effect of contraction on correlation tensor.—The analysis 
herein leads first to the changes in the spectrum tensor 
Tag*(k) TP ag?(x). Then the corresponding changes in the 
correlation tensor may be obtained from the Fourier trans- 
form relation, equation (6b): 


nos f Í f Tas etd e (E) (162) 
R.S()— f f f — (16b) 


In succeeding paragraphs I'44P(x) will be determined in terms 
of the initial spectrum tensor I'g*(k) for various cases. 

Spectrum tensors at upstream and downstream stations 
in terms of Q,.—In an earlier discussion, the Fourier coeff- 
cients Qae were chosen so as to define a field of turbulence 
confined to a large parallelepiped of volume 7, and vanishing 


everywhere outside; for this case lim ŠT Q.Qs* was to be 


identified with the correlation tensor Tag. For station A 
upstream of the contraction it will be convenient to special- 
ize this parallelepiped to a cube of edge D. Such a cube 
will, however, be distorted into a parallelepiped of edges 
LD, LD, D by the stream contraction by the time it reaches 
station B downstream. (See fig. 1(b).) The spectrum 
tensors for stations A and B, respectively, are therefore 
Dap*(k)= lim 5 TF - QAER E) 
(17) 


Tag Bn) = lim BET Qa" (x)Qs?*(x) 


Evaluation of spectrum tensor at downstream station, 
general case.—The identifications made in the last para- 
graph allow the spectrum tensor to be evaluated at station 
B in terms of the spectrum tensor at station A and the param- 
eters l, L, Ll defining the stream contraction between sta- 
tions A and B. For a single plane wave $Q,4— QA e'*z, 
which is transformed by the contraction into ĝa =Q eE, 
equations (13) give 


dtm. (OA) 1S) 


In the Fourier integral representation Ẹa* is to be interpreted 
as dga*, Ga” as dga”, Q^ as Q.“d7(k), and Q.? as Q.Pdr(x). 


Accordingly, 
Q= ae A Q,^ PRA Qa“ kgk a 


since lilds=dr(k)/dr(x). Thus 
QQ = L ka d "IE 49» (- ME ks 


Q,** OA buta a y Qr QI krkskakı (Qa A A* T. vksk a n | 
a 


ly l? x? 


(18) 


The corresponding relation kive the postcontraction 
and precontraction spectrum tensors is, by virtue of equa- 
tions (17), 


Dey! (E) kl (or i(k) kak 
Y 
T4344 (E) k, lese k 
D aaa : ratty (19) 
where k is related to x by 
ky, ke, k=l, kıy, loko, ls (20) 


Special case: axisymmetric contraction but arbitrary 
initial spectrum.—When the contraction is axisymmetric 
(L=), the equation of continuity in the form 


>k (21) 
Y 
may be used to simplify equation (19). The result may be 
written 
Id pa TAG) E ketTia ika) kı C 


ekt ks? -+k 


Tutk: kaks (1 — > | 
(ek +k? t ks y 


where the ratio «=1,*/l,7; for a large speed gain in the con- 
traction «<1. 

Special case: axisymmetric contraction, isotropic initial 
spectrum.—A further simplification occurs when the tur- 
bulence at station A is isotropic. In that case, the spectrum 
tensor T'as” (x) downstream of the contraction can be expressed 
explicitly to within an unknown multiplicative factor G(k). 
This results from the fact (ref. 3) that Ll'a41(£) must then be 
an isotropic second-order tensor; the isotropic property 
together with the continuity condition, equation (21), 
edes T'ag*(k) to be of the.form 


Tag (E) — G(E) (IG8.5—hL hg) (23) 
where 
OER | 1 for «=$ 
«8 "l0 for a*B 
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The right-hand side of equation (22) may be evaluated 
by means of equation (23). The diagonal terms reduce to 
relatively simple forms: 


(k?—k,*)k* 


-u 
PHOT, 00 CREE ES ia 
T227(x)=1, G(k) LEE CRM 

k kg (E — k 1—6) (25) 


(els? +k? +k) 
apd I'337(x) is obtained from Iss?(x) by replacing ks by k, and 


vice versa. The‘relation, equation 20), between k and x 


applies here. 

One-dimensional longitudinal spectrums.—If the form of 
the intitial spectrum tensor I'.4g4(£) is known, the correspond- 
ing one-dimensional spectrums (kj) can be calculated, 
according to the defining equation (7) as applied at the 
upstream station A: 


PAx2 f Í TA, ka, Eod Ead E, 


A particular case of isotropic turbulence is of special interest 
(ref. 4): In equation (23) for T.g4(k) the function G(k) is 
taken to be N(K4-437) ?, where Na are constants. "Then 


N(ks* + ks?) 
Ces a 


(26) 


and after integration 

_ aN 

Fr) q 345 
This one-dimensional longitudinal spectrüm is of the same 
form as an empirical relation obtained in reference 5 for tbat 
of isotropic turbulence in the initial period; this agreement 
is the special virtue of the form assumed for G(k). 

The one-dimensional lateral spectrum functions cor- 

responding to the same G(k) are readily evaluated; they 


are 
aN(3ki? Fy) 
rF 


The equality of the Fz and the F; functions results, of course, 
from the isotropy of the turbulence. 

The effect of the stream contraction on these one-dimen- 
sional spectrums is found by employing the postcontraction 
value of Taa, that 1s, the value Taa? appropriate to the 
downstream station B. Since Taa? is & function of the local 
wave-number vector x at station B, the equation correspond- 
ing to (26) is 


^ (6) =F )———7 3 rum 


P=? f [Peon K2, Ks). Kod Ks 


which is a function of x. 

For performing the integration and making later com- 
parisons of spectrums, it is convenient to transform from 
Ki; Ka, K3 to l5, ka, ls, where 

i — hl 
ka — Kal 
ky ks]l; 


and to define P, (k)-- L7 F,2, such that Í * Fd, = 
Í Fd; thus 


ky k: kg 
ne uem 


i TEN f VAT 


The spectrum tensor elements Faa” following an axisym- 
metric contraction have been evaluated in equations (24) 
and (25). With these values inserted and G(k) specified as 
before, the integrations of equation (25) are best effected in 
polar CODEOLUE las: The results are expressed most simply 
in terms of a ‘normalized’ longitudinal wave number k,/y 
as incorporated in the two parameters 


dksd lk; (27) 


2 
pee og 


The final result for the one-dimensional longitudinal spec- 
trum following an axisymmetric contraction (=l) 1s 


N 
FP (hi) =p 


| 24er 


(2+4e+2r+e+52) In (252) | (28) 


The corresponding result for the one-dimensional lateral 
spectrums following an axisymmetric contraction is 


zN ((8s—2)5 


P (E) FP (9 5 s (s 


10—96—0 PEE n +) 


ee Pn G) em 


and for «<1 (large speed gain) a simple but very close ap- 
proximation is | 


aN 1+2k, i? 
255 (FRA O O 


(The corresponding approximation for Fj?(k,) is not simple 
enough to warrant noting.) The parameters l, =la, and 
e=/,7/l;? are related to the initial and final Mach numbers of 
the main stream, if the fluid is air, by the equations 


Fy? (ky) = 


(az) GFS) 
Laz a) (30a) 


x.) SL) 


SPECTRUM OF TURBULENCE IN A CONTRACTING STREAM 107 


For incompressible flow (Mg, M4—0) these reduce to 
h= Usf Ua 


h="? 
e=1,73 


These postcontraction spectrums, equations (28) and 
(29), are compared with the initial spectrums in figures 2 
and 3, respectively; the comparison is based on an assumed 
initial stream Mach number of 0.05 (station A) followed by 
an axisymmetric contraction such that the final Mach num- 
ber is 2.0 (station B); the corresponding parameters are 
1,=29.8, L=0.382, e=0.00016. Consider first the longi- 
tudinal spectrums, figure 2. Normalizing factors are in- 
corporated such that the areas under the two curves, if re- 
plotted on a linear scale, would be the same; this normali- 
zation serves to differentiate changes in shape from changes 
in amplitude. The figure exhibits a rather striking distor- 
tion of the spectrum after traversing the stream contraction: 
the peak spectral density is shifted from zero wave number 
to kj/y—1.4 along with a general shift of density to the 
higher wave numbers. Associated with this change in shape 


is a reduction in amplitude by the factor ug?/u42, us? and uj? 
being the respective integrals of the spectral density curves. 
These integrals are evaluated in a later section. 


sia E EM 


Upstreom spectrum 


(30b) 
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ra 
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FIGURE 2.—Comparison of one-dimensional longitudinal spectrum upstream (F$) and 
downstream (FP) of axisymmetric contraction. Curves normalized to same area. Iso- 
tropic initial turbulence with G(k)—=N(k*++-73)7. Mi=0.05, M p-2.0, corresponding to 
120.821, «=0,0001637. 


One-dimensional longitudinal spectral densi 


The corresponding comparison for the lateral ona diim: 
sional spectrums is made in figure 3. In this case the axi- 
symmetric contraction has made very little distortion in the 
spectrum. ‘There is again a change in magnitude (this time 
an increase) in the ratio v5i/v,7. 

The changes in magnitude (that is, the changes in area 
under the spectral density curves) correspond to the changes 
un juat and v,?/v4? in the mean-square components of turbu- 
lence and are, at least qualitatively, well-known. The pre- 
dicted changes in the shape of the spectrum curves are ap- 
parently new. 

In the above comparisons both precontraction and post- 
contraction spectrums have been expressed in terms of the 
precontraction longitudinal wave number kı, whereas the 
local postcontraction wave number is =4,/l,. Consider, 
however, a representative longitudinal wave which has the 
form cos Iz at station A and cos mé at station B. If and £ 
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FIGURE 3.—Comparison of one-dimensional — spectrum upstream (F$) and down- 
stream (F7) of axisymmetric contraction. Curves normalized to same area. Isotropic 


initial turbulence with G(E) = N(k?--32)-*; M 490.05, M p2.0, corresponding to 7,229.821 
«e=0.0001637. 


are identified with the respective distances swept in time t 
by the moving waves over stationary hot-wire probes at 
stations A and B, respectively, then 4,e—k,U,t and r= 
Ky Ui li Uat. 
` 1 

hot wire in both cases is k,U,/27. The comparison based 
on k, therefore constitutes, in effect, & comparison of the 
time spectrums that would be seen by stationary hot-wire 
probes, in contradistinction to the space spectrums 
discussed in the earlier parts of the paper. ý 


Thus the (temporal) frequency seen by the 


EFFECT OF CONTRACTION ON MEAN-SQUARE VELOCITY COMPONENTS 
FOR ISOTROPIC TURBULENCE 


The mean-square velocity components of the turbulent 
field may be identified as the diagonal terms of the corre- 
lation tensor R.s(r) with r set equal to zero. Thus 


V-R,. (0) 
v= Ra(0) 


1 


where u, v, w have been written for qi, q», qs, respectively. 


. The evaluation of these means is much less laborious than 


the evaluation of the general correlation tensor. ln par- 
ticular, the evaluation of the ratio of the means Aug! fta, 
etc., may be made when the initial turbulence is specified 
to be isotropic- but no further details of its spectrum are 
known. These ratios will be calculated in the’ following 
paragraphs. 

Evaluations of u? and v! at upstream aaa 
to equation (6) 


Ey*(0)—u4 zz f f frs Atk) dr ON 


For isotropic turbulence I'zg has the form specified in equa- 
tion (23), whence 


us f f f 4h) hdr) 
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where G(k) is an arbitrary function. It is convenient to 
transform to spherical polar coordinates: 


b,—k cos Q. . 
k;—k sin @ cos 9 
(31) 
k.—k sin 0 sin p 
dr(k)=k? sin 0 dô dodk 
Then 
ug [ Gd) LL do | sin*o de (32) 


For the present purpose the function G(k), which, together 
with the condition of isotropy, defines the turbulence, may 
be left unspecified; the integral involving G(k) will cancel 
out in forming the ratio us?/u,?. Let this integral have the 
value H; then 


By virtue of the assumed isotropy 
wi. eH 


Evaluation of ratio of u? at downstream station to u? at 
upstream station.— The mean value up is obtained from an 
integration involving the spectrum tensor after the latter 
has been transformed by passage of the flow through the 
tunnel contraction; according to equation (16b) 


Ry? (0)=u?= f ff Ti? (x)d (x) 


For the present case, where the spectrum tensor at station 
A is assumed isotropic and the contraction is axisymmetric, 
the transformed tensor Lj?(x) has been determined in 


equation (24). Thus 
anit f f GHk (k’—k®dr (0) 
B ENS (ek? H- E tk) 


Because of the unspecified function G(k), it is convenient to 
change the variables of integration from the- components of 
x to the components of k. In other words, a transformation 
is made from the ‘‘wave-number space" of station B to the 
"wave-number space" of station A. The transformation 
. follows from the Cartesian relations 


dr (k)- d kid ksd E; 
dr (4) — d. ad. Kod xy 
_ ak, dk: dks 


hh h 


(33) 


_ isotropic. 


together with =l for an axisymmetric contraction, whence 


—, 1 ff PGK dr (E) 
“Ar JJ (e Ej? tkt kry? 


Again the polar-coordinate transformation (equation (31)) 
is made, with the result 


3 odo 
=A k*G (k)dk [^ ds r z a LA 


The first two integrals occur also in 14? (equation (32)), and 
they cancel in obtaining the ratio us us! fua; thus 


uy 3 3 JA sin? 0 dé 
ust 41, Jo (e cos? 0+sin? 6)* 
The final result may be written 


—] 2—e 
(1— ey 


us 3 


=> 4 
uUa? 4l? (3 ) 


tanh- T= | 


and an asymptotic expansion for small e is 


2 
=E= aa ae eT- (1 3- e)n. zF0(In 9| 
Us 4l, 

Equation (84) gives the ratio of the mean-square longi- 
tudinal velocity fluctuation downstream of an axisymmetric 
tunnel contraction to the corresponding mean square up- 
stream of the contraction, when the initial turbulence is 
The contraction is characterized by an increase 
in the stream speed in the specified ratio /, and a decrease 
in the lateral dimensions in the specified ratio l2; the param- 
eters l, L, and eeL?/l? are completely defined by the initial 
and final Mach numbers of the stream according to equa- 
tions (30a) and (80b). 


' 

The variation of u;/u,2 with the speed. ratio J; is plotted 
in figure 4 for two examples: in the first the flow is assumed 
compressible with & Mach number 0.05 at the start of the 
contraction; in the second the flow is assumed incompres- 
sible (M4, Ma—0). The Mach number scale at the bottom 
applies only to the compressible case, the J; scale to both 
cases. The salient characteristic of the curves is the marked 
reduction in the longitudinal component of turbulence with 
increasing speed ratio l. 

Compressibility is seen to have but a secondary effect, 
which is appreciable only at supersonic speeds. Note 
(equations (30) and (34)) that with /, as the independent 
variable, the effect of compressibility appears only in the 
parameter e. The physical significance of e follows from tho. 
definition of l as the speed ratio provided by the contraction 
and L? as the area ratio of the contraction (in the axisym- 
metric case considered), with e=}?/L?. For supersonic final 
speeds it is more proper to speak of a converging-diverging 
nozzle than a contraction, the term ‘contraction’ having 
been retained herein primarily for reasons of past usago. 
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The basis of the compressibility effect may be summed up 


in the following way. The influence of an axisymmetric 
stream contraction arises from distortion of the fluid ele- 


ments, as described by the parameters /, and 4L. (See fig. 1 (b).) 
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FIOURE 4,—T ypical examples of selective effect of axisymmetric stream contraction without 
decay on components of turbulent Intensity, showing influence of compressibility. Iso- 
tropie turbulence at l1. 


These parameters are related by the continuity condition 
oh ly) = 1, where o is the density ratio. Thus compressibility, 


in allowing c to deviate from unity, changes the relation ` 


between l and /; somewhat, and consequently modifies the 
contraction effect. 

The graph of equation (34) in figure 4 is primarily for 
illustrative purposes; a form more useful for engineering 
applications is given in figure 5. The single curve provides 


the variation of Vus fun? with both /, and e; lh and « may 
be determined from the initial and final Mach numbers by 
means of the simple relations (30a). 

Evaluation of ratio of v at downstream station to p! at 
upstream station.— The value of vs results from an inte- 
gration involving the transformed spectrum tensor, accord- 
ing to equation (16b) 


(0) == M f Í TP (9d «(9 


For isotropic initial turbulence and an axisymmetric con- 
traction, the transformed spectrum tensor T'x;?(x) has been 
evaluated ; in equation (25). ‘Thus 


2k? (1—e 


m= f f frof- MÀ erppr 
MES 
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FIGURE 5.—Varintion of relative root-mean-square longitudinal velocity fluctuation with 
both speed ratio Jı and area ratio /2 for axisymmetric contraction and isotropfe turbulence 
at i41. 


Again it is convenient to transform from x-space to k-space 
(equations (33)) and to introduce polar coordinates k, e, 
and 0 (equations (31)). ‘The integrations with respect to k 


and e are readily disposed of, with the result 


: * T ain? 2 
opal 2 f sin? 6d8—2x(1—«) f Di dcos Cdi y 


sin? 8+ e cos? 8 


x(1—e)? sin’ 0 cos? 8d0 ] 


o (sin? 0- e cos? 6)" 
where H= ji i kiG(k)dk, as before. Upon carrying out the 


integration and dividing by v EE 3 TH, there is obtained 
finally 


]—e (39) 


He tanh"! yi= | 


For small e this has the asymptotic expansion 


ve 


s Sil tery th pero 


Equation (35) gives the ratio of the mean-square lateral 
velocity fluctuation downstream of an axisymmetric tunnel 
contraction to the corresponding mean square upstream of 
the contraction, where the turbulence has been assumed to 


be isotropic. The variation of 4v KT with the speed ratio 
l, is plotted in figure 4, which already contains the graph of 


Jusiju d discussed earlier; again the two cases are incom- 
pressible flow and compressible flow with an initial Mach 
number of 0.05. For /,21 and incompressible flow, the ` 
lateral component of turbulence is seen to increase steadily 
with l}, in marked contrast to the decrease exhibited by the 


Pad 
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longitudinal component. The curve (of the lateral com- 
ponent) for compressible flow begins to differ sensibly from 
the curve for incompressible flow for downstream Mach 
numbers above 0.3; above sonic speed compressibility is seen 
to effect a complete reversal of the curve. The over-all 
effect of compressibility on the contraction effect is thus 
much greater for the lateral than for the longitudinal com- 
ponent of the turbulence. 

The graph of equation (35) in figure 4 is primarily illus- 
trative; a form more useful for engineering applications is 
given in figure 6. The single curve provides the variation 
of Jos va with both l and e; l and e may be determined 
from the initial and final Mach numbers by means of equa- 
tions (30a). 

- DECAY CONSIDERATIONS 


CRITERIA FOR NEGLIGIBLY SMALL DECAY 


The basis of the present analysis of the contraction effect 
is embodied in equations (118) relating the precontraction 
and postcontraction vorticity distributions. The simplicity 
of this result and its derivation arises from the neglect of the 
turbulent decay; by decay is meant the viscous dissipation 
and all the (nonlinear) intermixing processes of the eddies 
which together cause the mean turbulent intensity to dimin- 
ish with time. The postulation of an inviscid fluid elimi- 
nated the viscous dissipation, and the limitation to very 
weak turbulence eliminated the intermixing processes. 
(While there can be no dissipation in an inviscid fluid, the 
intermixing processes ordinarily associated with decay will 
occur.) In order to assess the influence of these assumptions, 
equations (11a) will now be derived in a more general fashion 
with the Navier-Stokes equations as the starting point. For 
simplicity the fluid is taken to be incompressible, since the 
major conclusions are unaffected thereby. 


General formulation of changes in vorticity. By rearrange- | 


ment and cross differentiation to eliminate the pressure term 
(ref. 8, p. 578), the Navier-Stokes equations can be trans- 
formed into ‘ 


Do oq 


Di” dz, (36) 


Ò 
+ we 2n + ws an Fa 
T2 


and two similar Squatións, where w=w;, «vs, w; 18 the vor- 
ticity and g’=q;’, qs’, gs’ is the resultant velocity. Now let 
q' be the sum of a stream velocity U, V, W and a turbulent 
velocity field q=qı, qa, q3; also, let curl U, V, W=0, so that 
e is just curl q. Then equation (36) becomes (n iensor 
notation) E 
D mE VET 3z T8 ID Eta 
Ru uc ANE 


Contraction Decay 


(37) 


and there are two similar equations. The first set of terms 
on the right-hand side is identified as the contraction effect, 
the second set as the decay effect. First the decay terms 
will be neglected in an attempt to recover equations (11a); 
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pie 6— Variation of relative root-mean-square lateral velocity fluctuation with both 
speed ratio l; and area ratio {2 for axisymmetric contraction and isotropio turbulenco at 
lhal, 


then the neglected decay terms will be examined and cri- 
teria for their neglect arrived at. 

Neglect of decay terms.—Equation (37) minus the viscous 
terms reads, in expanded form, 


Do _ AL) o(U-rgq), -DU +a) 
De Q1 H+ ey O23 +s Og (38) 


In this and the earlier equations $5 is the 'Lagrangian? de- 


rivative following the fluid motion. Now consider a line 
segment à, 522, 523 following the fluid motion: its Lagrangian 
derivative can be shown to be 


a br X0) us, XD s, Wta (39) 


and two similar equations. It can be seen that a solution of 
equations (38) and (39), together with nid companion 
SHURDONS is given by 


Q1, 0, Q3 77 02, 023, bX 


for all time 7; this result is well known. Now complete the 
neglect of the decay terms by omitting the terms in q, in 
equation (38) and correspondingly in equation (39). By 
this neglect the turbulent perturbations of the flow stream- 
lines have been suppressed: this can be inferred from the 
revised equation (39). If the particles are at station A at a 
time t=0 and reach station B at time t=t, there results 


. . . ôx? « . 62:4? ‘ 
and two similar equations. But zaa 18 Just li, $z,4 8 ls, 
s 1 2 
d 9 is 1, 
Ts 


Therefore equations (lla) have been recovered for the in- 
compressible case (density ratio c=1). 

Consideration of inertial decay terms.—In equations (37) 
to (39) the decay terms not involving » are the inertial or 
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intermixing terms. These are seen to be nonlinear. The 
condition for their neglect is evidently l 
Ogi oU 
Drs <j weg ozs (40) 


and two similar conditions between q and V, q, and W, 
respectively. In & contraction like that of a wind tunnel, 
the dominant velocity gradients will be te? or and 
3 
these wil be of the same order of dE magnitude. A 
sufficient condition to replace expression (40) is therefore 


bU 
Ox; 


Oda 
OZg 


(41) 


that is, all of the turbulent velocity gradients are very much 
less than the axial gradient of the stream velocity. This is 
essentially the assumption underlying the distortion equa- 
tions (10), which led directly to the vorticity changes (11a) 
in Taylor's method. 

In statistical terms an approximate inference Pon equa- 
tion (41) is, for isotropic turbulence, 


Y ($=) 


But by definition of à this may be written 


Wi oU 
<E 


The “microscale” à may be interpreted as a sort of average 
eddy diameter weighted in favor of the smaller eddies. 
Equation (42) may be accepted as a practical criterion. for 
the neglect of the inertial decay terms, equivalent to one of 
the two assumptions underlying equation (10). The other 
assumption, neglect of viscosity, is considered next. 

Consideration of viscous decay term.—The viscous decay 
term in equation (37) is the term containing v. This term is 
linear and so will affect individual plane waves separately 
without mutual interference. The magnitude of the term 
may be estimated to a sufficient approximation by considering 
a wave carried along by the contracting stream and neglect- 
ing (for this term only) the distortion of the wave imposed 
by the contraction. Thus & component of the wave may be 
written 


«52 


(42) 


w= Qie! (ez 51UD 


Then, if the inertial decay terms of equation (37) are negli- 


oU| |oU| /||oU 
gible, the equation reads, with oz, azl . oz: 
E On, — Vui 
E “i LRA E AE E 


Accordingly, viscosity may be neglected for that portion of 
the spectrum which satisfies the inequality 


rhe (43) 


In the "initial period” of decay, if the inertial decay criterion 
(42) is satisfied, the major part of the speetrum will satisfy 
condition (43). 


ROUGH ESTIMATION OF MUTUAL EFFECTS OF DECAY AND CONTRACTION 


When the decay effects are not negligible compared with 
the contraction effects (see criterion developed in the last 
section) the theoretical basis of the present theory of the 
contraction effect is violated. Because negligible decay is 
more the exception than the rule, there is considerable 
incentive to attempt to apply the theory outside the valid 
range by means of assumptions concerning the simultaneous 
effects of decay and contraction. 

Suppose, now, the decay and the contraction are considered 
to occur alternatively.in small steps, starting from isotropic 
turbulence. Each stream tube is considered to contract 
stepwise: between steps there is decay without contraction; 
at each step there is a sudden contraction without decay. 
Let the change in speed ratio per step be dl, the reduction 


in v? due to decay be (dup, and the reduction in v? due to 


contraction be (duje. Express the effect of decay in the 
absence of contraction in the form 


(gos 


where l is a function of the time of travel (decay time) t, 
and the effect of contraction in the absence of decay in the 


(44) 


form 
u? 
(4) =c () . i (45) 
The corresponding differential forms are 
(du?)p DL) 
Ec Xx Dh) dh (46) 
(df) je. C” r, dl, (47) 


a Oh) 


The assumption is now made that equation (46) applies to 
the decay effect per step and equation (47) to the contraction 
effect per step, the only interaction being in the common 
ui' This assumption neglects the tendency of the decay 
process to counteract the anisotropy ‘introduced. by the 


4 It is known that In the ‘initial’ period of decay T, — ud. Equation (46) amounts 
K S 7 


to replacing the —u? on the right-hand side by —(t?) desay oniy; gome defense may bo 
made of this approximation, considering the progressive deviation :rom isotropy. 
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contraction effect. The total effect per step is then 


dé ford), D'Q) | 
LOW +a) | aa 


whence upon integration the over-all effect 1s 


+ —3 


=D (48) 


That is, if the effect of contraction alone is expressed by 
C(l) (equation (45)) and the effect of decay alone by D(h) 
(equation (44)), then the joint effect under the assumption 
is expressed by the product C(L) D(L). 

Equation (48) is intended to provide a very rough adjust- 
ment of the theoretical contraction effect C(l) to account 
for decay. This adjustment will be made in the attempt 
to compare the theory with experimental results in which 
the decay effects are of the same order as the contraction 
effects. 

Equation (48) refers to the longitudinal velocity com- 
ponent v; an equation of the same form is obtained for the 
lateral component ve. For both cases the function D(L) is 
taken to be the right-hand side of the empirical decay law 


(ref. 9) 
(=) -——=— 
Ua?) p Yu 


140.58+— t(1) 


= D() (49) 


for isotropic turbulence in the initial period. The decay 
time ¿(l) in the formula is the time required by'& particle 
of the main stream to pass through the contraction, the 
initial velocity being Ua and the final velocity l Ua. 


COMPARISON WITH EXPERIMENT 


There appear to have been no experimental investigations 
with which to compare the predicted changes imposed by a 
stream contraction on the spectrum of the turbulence, or 
‘on the correlation tensor of the turbulence. The available 
experimental data seem to be limited to measurements 


bearing on the changes in the root-mean-square velocity | 


components. These data apply, moreover, to conditions 
outside the proper scope of the present theory in that large 
decay effects are present. "The experimental data are'there- 
fore compared with a crude extension of the theory in which 
the decay is allowed for in first approximation. (See pre- 
ceding section.) 

The most extensive data are those of MacPhail (ref. 10), 
which in effect cover a range of contraction ratios from 
l1 to /,29.65 inasmuch as measurements were made at 
various stations along the contraction. Isolated points for 
particular contraction ratios were obtained from investiga- 
tions made for other purposes by Dryden and Schubauer 
(ref. 11) and by Hall (ref. 12). Only those points were 
chosen for which the initial turbulence was indicated to be 
approximately isotropic. 
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In table I are listed, for the three experimental arrange- 
ments, the parameters used in the estimation of the decay 
factor (equations (48) and 49)). In reference 10 the initial 


‘stream velocity U, and scale of turbulence La were given. 


In reference 11 the value of U, was given, and the value of 
L, was taken to be 0.05 foot, the only scale mentioned; it 
was not clear, however, whether this value of scale applied 
with or without screens. In reference 12 the value of U, 
was inferred from collateral information and is somewhat 
uncertain; the scale La was estimated from the dimensions 
of the honeycomb. In all three experimental arrangements 
the initial relative levels of turbulence were specified. The 
decay time ¢ of the turbulence was computed as the time for 
a particle to traverse the contraction; the value arrived at 
for Hall’s data (ref. 12) reflects the uncertainty in the 
assumed U,. 


Root-mean-square longitudinal velocity components.— 
The comparison of the theory, including estimated decay, 
with experiment for the longitudinal component of turbulence 
is given in figure 7. The theoretical curve, in each instance, 
is the product of a value computed for contraction along, 
neglecting decay, (obtainable from fig. 5) and a value esti- 
mated for decay alone neglecting contraction. (See equa- 
tions (48) and (49).) The agreement with MacPhail’s data 
and with Hall’s single point can be considered good. The 
agreement with the Dryden-Schubauer point, on the other 
hand, is poor; a slight.improvement would result on correc- 
tion for the spurious contribution of the noise background. 


Experiment Pie el 


o MacPhail (ref. 10) 
© Dryden-Schubouer (ref. 11) 
o Holl (ret 12) 


Theory with estimated decoy 
for experiment of MocPhoil 
—-— for experiment of Dryden-Schubauer 
for experiment of Holl 


Relotive rms longitudinal velocity fluctuation 
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FIGURE 7.—Comparison of predicted axisymmetric contraction effect with oxperimont for 
longitudinal component of turbulence, with decay allowed for In first approximation. 
Initial isotropic turbulence assumed. 
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Root-mean-square lateral velocity components.—Com- 
parison of the theory, again including estimated decay, with 
experiment for the lateral component of turbulence is given 
in figure 8. There is complete disagreement with MacPhail’s 
data and Hall’s single point, and on the other hand, good 
agreement with the Dryden-Schubauer single point. Thus 
there is the curious result that MacPhail’s and Hall’s data 
agree well with theory for the longitudinal component and 
disagree entirely for the lateral component, whereas the con- 
verse is true for the Dryden-Schubauer data. 

Discussion.—The uncertainty both in the manner of 
estimating the decay effect and in the data (table I) on which 
the estimate was based is still far from sufficient to account 
for the discrepancies between theory and experiment for the 
lateral component of turbulence. The very large amplifica- 
tion found by MacPhail is particularly hard to explain. On 
the other hand, the experimental data of the several ob- 
servers show considerable disagreement, especially when 
differences in decay are allowed for. This disagreement 
would tend to cast doubt on the validity of some of the data; 


o MacPhail (ref. 10) 
© Dryden-Schubouer (ref. 11) 
a Hall (ref. 12) 


Theory with estimated decay 
for experiment of MacPhail 
—-— for experiment of Dryden-Schubauer 


Relative rms loteral velocity fluctuation 
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Fiaure 8,—Comparison of predicted axisymmetric contraction effect with experiment for 
lateral component of turbulence, with decay allowed for in first approximation. Initial 
fsotrople turbulence assumed. 


the disagreement may also be in part a consequence, pre- 
dicted by the theory, of possible differences of the initial 
spectrums from each other and from isotropy. 


CONCLUDING REMARKS 


The original aim of this paper was to provide a quantitative 
explanation of the observed changes in ili: root-mean-square 
velocity components of the turbulence of & wind-tunnel 
stream after passing through the tunnel contraction. The 
simplifying assumption of negligible decay was made to make 
the analysis tractable, although the decay and contraction 
effects are ordinarily comparable. "The analysis on this basis © 
disclosed, in addition to the above integrated effects, pro- 
nounced changes in the spectrum of the turbulence. The 
changes in the shape of the spectral density curves, as dis- 
tinguished from over-all changes in amplitude, would appear 
to be considerably less sensitive to modification by decay than 
would the mean-square velocity components. For this 
reason, and because such spectral changes have not pre- 
viously been discussed, the emphasis of the present paper 
has been placed most heavily on these spectral effects. 

In particular, it has been found that the one-dimensional 
longitudinal spectrum for isotropic turbulence exhibits a 
rather interesting change in shape downstream of the con- 
traction; the center of gravity of the curve of spectral density 
versus longitudinal wave number is shifted substantially to 
higher wave numbers, the resulting distortion moving the 
peak of the curve well to the right of its initia] position above 
the origin. The distortion is quite pronounced and would 
appear to be readily amenable to experimental observation. 

The restrictive assumption of negligible decay largely de- 
feats the origina] aim of the paper. Nevertheless, for 
practical reasons an attempt has been made to provide a 
crude extension to the theory in which decay is allowed for in 
first approximation. With this approximation the theory 
has been compared with experimental values of the contrac- 
tion effect on the longitudinal and lateral component root- 
mean-square velocity fluctuations. The agreement for the 
longitudinal component is good, whereas there appears to 
be almost complete disagreement for the lateral component, 
the experimental data themselves being in conflict. It is 
perhaps premature to attempt any general conclusion. For 
the present, the theory as augmented by the estimated decay 
effect may be useful in wind-tunnel-design applications. 

It is clear that the tunnel contraction effect on the com- 
ponents of turbulent intensity cannot be represented by fixed 
fractional changes independent of the character of the initial 
turbulence. Instead the separate factors for the longi- 
tudinal and lateral components depend markedly on the 
spectrum of the turbulence. For initial isotropy, however, 
unique factors are predicted that, when decay is neglected, 
are independent of the details of the spectrum. — 


Lewis Fricut PROPULSION LABORATORY 
NATIONAL ÅDVISORY COMMITTEE FOR AERONAUTICS 
CLEVELAND, Onto, August 30, 1951 
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su GF RI APPENDIX A 
! | ' SYMBOLS 
The following notation is used in this report: ; ; ( (e—1) £j? ) 
tion (28) ( === _ 4 
The subscripts 1, 2, 3 refer to a rectangular coordinate : poremouny in equation (28) y 
system with the 1-axis alined with the axis of the main flow | / time l 
and directed downstream, the 2-axis directed horizontally, U main-stream velocity l 
and the 3-axis directed vertically. Separate systems are used | %?;Y= 41, d» ds disturbance velocity components 
with the origins at stations A and B, respectively. (See fig. 1.) X length of wind-tunnel contraction (distance be- 


Vector and tensor notations are used interchangeably; for 
example, k=a= (kı, ke, k3), where a=1, 2, or 3 designates 
a vector with components £E;, ke, and ks. 


C(l) function defined in equation (45) 

DU) function defined in equation (44) 

D edge length of cube within which turbulent field 
is defined. 

e base of natural logarithms 

r= In, F>, or F; 

F, one- -dimensional longitudinal spectral oss 
(see equation (7))' 

Fa, F, one-dimensional lateral spectral densities (see 
equation (7)) : 

G(k) function appearing in isotropic spectrum tensor 

Ws. — oganstunt ( i ka (k) dk) 

Im imaginary part of i 


= 4—1 
, Ks, Byki thy’, ky, katk’, respectively 


K, 2 

k amplitude of k(=/k?+k2+k,*) 

k= ka= (ki, ke, k3) wave number vector (station A) 

L lateral scale 

la li, b, A 

l stream velocity at station B divided by stream 
velocity at station A (see fig. 1) 

ls stream breadth at station B divided by stream 
breadth at station A (see fig. 1) 

ls stream height at station B divided by stream 
height at station A (see fig. 1) 

M Mach number of main stream 

amplitude of special isotropic spectrum tensor 

(see following equation (26)) 


Q= Qa= (Q1, Qs, Q9) disturbance wave amplitude vec- 
tor 

q= = (qı; d», qs) disturbance velocity vector 

Raet) seilicet tensor (ref. 6 ) 

Re real part of 

r magnitude of r=-+/r?-+7.?7-+7;? 

r— Ta™= (ri, rs, 73) Separation vector of two correlated 
points 


69 


2 
parameter in equation (28) (=+ 1) 


a aM 


tween stations A and B) 

used occasionally in place of zı 

La= (zi, 25, x3) position vector (station A) 

aplk) spectrum tensor (ref. 3) 

constant in special isotropic spectrum tensor 
(see following equation (26)) (y=1/Z) 

contraction parameter (=/,7/l,?; see fig. 1) 

alternating tensor defined after equation (9) 

polar angle (equation (31)) 

magnitude of x( 9x? + i Es?) 

Ka™ (Ki, K2, K3) transformed wave number vector 

` (station B); (xa=kafla) 

kinematic viscosity 

magnitude of (=V Fi F4) 

a= ($n &, &) transformed position vector 
(station B) (see equation (10) and fig. 1(b)) 

summation over a for a=1, 2, 3 


Ix * D ^ ~ HI NH 
i g 


|vrr re & 


stream density at station B divided by stream 
density at station A 

a volume 

azimuth angle (equation (31)) 

Q2, — (Qi, R, Q3) vector amplitude of vorticity wave 


Wa= (w, w, w) vorticity vector 


IE T 1 
| ‘ 


Superscripts: 

A measured in vicinity of station A, upstream of 
contraction 

B measured in vicinity.of station B, downstream 
of contraction 

ii complex conjugate 


Subscripts: 
A measured in vicinity of station A, upstream of 
l contraction 
B measured in vicinity of station B, downstream 
of contraction 
. take on values 1, 2, or 3‘and designate tensor 
quantities 
specific values of a, B, y, or 6 


a,B,7,6 
1,2,3 


A symbol with the mark ~ above it refers to a single plane 
wave. <A bar over a symbol designates an average (usually 
a spatial average); a bar under a symbol designates a vector. 


SPECTRUM OF TURBULENCE IN A CONTRACTING STREAM 115 


APPENDIX B 
EQUIVALENCE OF SPACE AND TIME AVERAGES IN STATISTICALLY STEADY, HOMOGENEOUS TURBULENCE 


The definitions of statistical homogeneity and statistical 
time independence will first be made precise. Let F' (z, y, z, t) 
be some property of & turbulent field that varies in time 
and from point to point; thus F may be the pressure, or any 
of the velocity components, or & correlation of velocity com- 
ponents at two points of fixed separation, (x, y, z) being one 
of the two points. If, for all choices of the property F, (a) 
the average of F over a time T— œ is independent of (z, y, 2), 
the turbulence is defined to be statistically homogeneous; 
if (b) the average of F over a volume V— œ is independent 
of t, the turbulence is defined, in the sense used herein, to be 
statistically steady or time-independent. The respective 
averages are supposed to be approached uniformly, in the 
mathematical sense, as 7' or V, respectively, approaches 
infinity. (A statistically steady or “stationary” condition is 
defined differently in the theory of random processes.) 

It will now be proved that if the turbulence satisfies the 
two conditions (a) and (b), the time and space averages 
defined therein are equal. In this proof, no resort will be 
made to the ‘ergodic hypothesis" of statistical mechanics, 
which leads to the equivalence of the time average and the 
“ensemble” average. The possibility of the joint existence of 
the conditions (a) and (b) probably amounts, however, to 
just as fundamental an assumption. 

The space average will be made over a parallelepiped of 
edges a, b, c and the time average over a time T, and then a 
limiting process will be applied. The average of F over both 
space and time is thus 


— : 1 T (^c Pb fa 
P= dm oop f Í f f Fdzdydzdt — (B1) 


Any order of integration is permissible, since the integra- 
tion limits are constants. If the time integration is per- 
formed first the expression may be written 


es . 1 c fb fa T 
Fa dim abe? Í Í Í ai Fdt ) dzdydz 


By virtue of the postulated uniform convergence of the time 
and space averages, the operation 7 lim ! may be brought under 
the integral sign: 


- T . 
T,= dim s IN i + (lim 7 E f Fat) dedydz 


= lim - |^ IN F,dzdydz 


where F, is the time average of F. But, by condition (a), 
F, is independent of z, y, and z. Therefore 


(B2) 


zl 


F= 


(B3) 


Alternatively, the space integration and limiting process 
may be performed first: 


Pa, P, lim an, (im, ais Sedo I Fdadydz) i 
af Fat 
= im 


where F, is the space average of F. By condition {b), F, is 
independent of t; therefore 


(B4) 


F, =F, (B5) 
Equations (B3) and (B5) together state that 
F=F,=F,, (B6) 


or the space average, the time average, and the space-time 
average are all equal. 
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